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Abstract 

In  this  paper,  we  show  the  relationship  between  two  algorithms  and  optimization  problems 
that  are  the  subject  of  recent  attention  in  the  networking  and  control  literature.  First,  we 
obtain  some  results  on  averaging  algorithms  over  acyclic  digraphs  with  fixed  and  controlled- 
switching  topology.  Second,  we  discuss  continuous  and  discrete  coverage  control  laws.  Further, 
we  show  how  discrete  coverage  control  laws  can  be  cast  as  averaging  algorithms  defined  over  an 
appropriate  graph  that  we  term  the  discrete  Voronoi  graph. 


1  Introduction 

Consensus  and  coverage  control  are  two  distinct  problems  within  the  recent  literature  on  multiagent 
coordination  and  cooperative  robotics.  Roughly  speaking,  the  objective  of  the  consensus  problem  is 
to  analyze  and  design  scalable  distributed  control  laws  to  drive  the  groups  of  agents  to  agree  upon 
certain  quantities  of  interest.  On  the  other  hand,  the  objective  of  the  coverage  control  problem  is 
to  deploy  the  agents  to  get  optimal  sensing  performance  of  an  environment  of  interest. 

In  the  literature,  many  researchers  have  used  averaging  algorithms  to  solve  consensus  problems. 
The  spirit  of  averaging  algorithms  is  to  let  the  state  of  each  agent  evolve  according  to  the  (weighted) 
average  of  the  state  of  its  neighbors.  Averaging  algorithms  has  been  studied  both  in  continuous 
time  [2,  3,  4,  5]  and  in  discrete  time  [6,  5,  7,  8,  9,  10,  11,  12].  In  [2],  averaging  algorithms  are 
investigated  via  graph  Laplacians  [13]  under  a  variety  of  assumptions,  including  fixed  and  switching 
communication  topologies,  time  delays,  and  directed  and  undirected  information  flow.  In  [3],  a  series 
of  consensus  protocols  are  presented,  based  on  the  regular  averaging  algorithms,  to  drive  the  agents 
to  agree  upon  the  value  of  the  power  mean.  A  theoretical  explanation  for  the  consensus  behavior 
of  the  Vicsek  model  [14]  is  provided  in  [7],  see  also  the  early  work  in  [6],  while  [5]  extends  the 
results  of  [7]  to  the  case  of  directed  topology  for  both  continuous  and  discrete  update  schemes.  The 
work  [8]  adopts  a  set-valued  Lyapunov  approach  to  analyze  the  convergence  properties  of  averaging 
algorithms,  which  is  generalized  in  [9]  to  the  case  of  time  delays.  Asynchronous  averaging  algorithms 
are  studied  in  [10].  The  work  [11]  analyzes  the  averaging  algorithms  in  the  framework  of  partial 
difference  equations  over  graphs  [15].  The  works  [16,  17]  survey  the  results  available  for  consensus 
problems  using  averaging  algorithms.  In  the  scenario  of  coverage  control,  [18]  proposes  gradient 
descent  algorithms  for  optimal  coverage,  and  [19]  presents  coverage  control  algorithms  for  groups 
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of  mobile  sensors  with  limited-range  interactions.  Also,  we  want  to  point  out  that  a  special  kind  of 
directed  graphs,  namely  acyclic  digraphs,  are  presented  in  the  literature  to  describe  the  interactions 
of  agents  in  leader-following  formation  problems,  e.g.,  [20,  21,  22]. 

The  contributions  of  this  paper  are  (i)  the  investigation  of  the  properties  of  averaging  algorithms 
over  acyclic  digraphs  with  fixed  and  controlled-switching  topologies,  and  (ii)  the  establishment  of 
the  connection  between  discrete  coverage  problems  and  averaging  algorithms  over  acyclic  digraphs. 
Regarding  (i),  our  first  contribution  is  a  novel  matrix  representation  of  the  disagreement  function 
associated  with  a  directed  graph.  Secondly,  we  prove  that  averaging  over  an  fixed  acyclic  graph 
drives  the  agents  to  an  equilibrium  determined  by  the  so-called  “sinks”  of  the  graph.  Finally,  we 
show  that  averaging  over  controlled-switching  acyclic  digraphs  also  makes  the  agents  converge  to 
the  set  of  equilibria  under  suitable  state-dependent  switching  signals.  Regarding  (ii),  we  present 
multicenter  locational  optimization  functions  in  continuous  and  discrete  settings,  and  discuss  dis¬ 
tributed  coverage  control  algorithms  that  optimize  them.  We  discuss  how  consistent  discretizations 
of  continuous  coverage  problems  yield  discrete  coverage  problems.  Finally,  we  show  how  discrete 
coverage  control  laws  over  the  discrete  Voronoi  graph  can  be  casted  and  analyzed  as  averaging 
algorithms  over  a  set  of  controlled-switching  acyclic  digraphs.  Various  simulations  illustrate  the 
results. 

The  paper  is  organized  as  follows.  Section  2  introduces  our  novel  matrix  representation  of  the 
disagreement  function,  and  then  reviews  the  current  results  on  consensus  problems.  We  also  present 
convergence  results  of  averaging  algorithms  over  acyclic  digraphs  with  both  fixed  and  controlled- 
switching  topologies.  Section  3  presents  locational  optimization  functions  in  both  continuous  and 
discrete  settings,  and  then  discusses  appropriate  coverage  control  laws.  The  main  result  of  the  paper 
shows  the  relationship  between  averaging  over  switching  acyclic  digraphs  and  discrete  coverage. 
Various  simulations  illustrate  this  result,  and  show  the  consistent  parallelism  between  the  continuous 
and  the  discrete  settings.  Finally,  we  gather  our  conclusions  in  Section  4.  For  easy  reference,  we 
review  some  basic  facts  and  standard  notations  from  nonsmooth  analysis. 

Notation 

We  let  N,  M+  and  M+  denote,  respectively,  the  set  of  natural  numbers,  the  set  of  positive  reals, 

O 

and  the  set  of  non-negative  reals.  For  any  set  5  C  l2,  we  let  S,  and  dS  denote,  respectively,  the 
interior,  and  boundary  of  S.  The  quadratic  form  associated  with  a  symmetric  matrix  B  £  Rnxn 
is  the  function  defined  by  x  e- >  xT Bx.  The  map  f  :  X  —>  Y  and  the  set-valued  map  /  :  X  z4  Y 
associate  to  a  point  in  X  a  point  in  Y  and  a  subset  of  Y,  respectively.  The  sum  of  m  subsets 
St ,  i  £  {1, . . . ,  m,}  in  a  vector  space,  denoted  by  Y1T=i  consists  of  all  vectors  of  the  form  YliLi 
where  st  £  Si,  for  i  £  {1, ... ,  m}. 


2  Averaging  algorithms  over  digraphs 

2.1  Preliminaries  on  digraphs  and  disagreement  functions 

A  weighted  directed  graph,  in  short  digraph,  Q  =  ( U ,  £,  A)  of  order  n  consists  of  a  vertex  set  U 
with  n  elements,  an  edge  set  £  £  2WxW  (recall  that  2U  is  the  collection  of  subsets  of  U),  and  a 
weighted  adjacency  matrix  A  with  nonnegative  entries  aij,  i,  j  £  {1, . . .  ,n}.  For  simplicity,  we  take 
U  =  {1, ,  ?r}.  For  i,j  £  {1, . . . ,  n},  the  entry  atj  is  positive  if  and  only  if  the  pair  ( i ,  j )  is  an  edge 
of  Q,  i.e.,  >  0  (i,j)  €  £■  We  also  assume  an  =  0  for  all  i  £  {1, . . .  ,n}  and  a^  =  0  if  (i,j)  0  £, 

for  all  i,j  £  {1, . . .  ,n}  and  i  A  j •  When  convenient,  we  will  refer  to  the  adjacency  matrix  of  Q  by 

AG)- 

Let  us  now  review  some  basic  connectivity  notions  for  digraphs.  A  directed  path  in  a  digraph  is 
an  ordered  sequence  of  vertices  such  that  any  two  consecutive  vertices  in  the  sequence  are  an  edge 
of  the  digraph.  A  cycle  is  a  non-trivial  directed  path  that  starts  and  ends  at  the  same  vertex.  A 
digraph  is  acyclic  if  it  contains  no  directed  cycles.  A  node  of  a  digraph  is  globally  reachable  if  it 
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can  be  reached  from  any  other  node  by  traversing  a  directed  path.  A  digraph  is  strongly  connected 
if  every  node  is  globally  reachable. 

Remark  2.1.  The  previous  definition  of  adjacency  matrix  follows  the  convention  adopted  in  [2], 
where  aij  >  0  <S>  (i,j)  £  £■  On  the  other  hand,  in  [16],  aij  >  0  (j,  i)  G  £■  This  difference  arises 

from  a  different  meaning  of  the  direction  of  an  edge.  In  [2],  a  directed  edge  ( i,j )  G  £  means  node  i 
can  ’see’  node  j,  i.e,,  node  i  can  obtain,  in  some  way,  information  from  node  j .  We  refer  to  this 
as  the  communication  interpretation.  In  [16],  a  directed  edge  (■ i,j )  G  £  means  that  the  information 
of  node  i  can  flow  to  node  j.  We  refer  to  this  as  the  sensing  interpretation.  The  difference  leads 
to  different  statements  of  various  results.  For  example,  having  a  globally  reachable  node  in  the 
communication  interpretation  is  equivalent  to  having  a  spanning  tree  in  the  sensing  interpretation. 


The  out-degree  and  the  in-degree  of  node  i  are  defined  by,  respectively, 

n  n 

don)  (i)  —  '  ajj ,  d\u  (i)  —  ^  Oj , . 

j= 1  i=1 

The  out-degree  matrix  Dout(Q)  and  the  in-degree  matrix  Dm{Q)  are  the  diagonal  matrices  defined  by 
(Ajut (£?))»»  =  G?out(*)  and  (Din(G))u  =  d-m{i ),  respectively.  The  digraph  Q  is  balanced  if  Dou t(Q)  = 
Dm{Q).  The  graph  Laplacian  of  the  digraph  Q  is 

L(Q)  =  Dout(G)  -  A(G), 


or,  in  components, 


3  =  h 

3  +  i- 

Next,  we  define  reverse  and  mirror  digraphs.  Let  £  be  the  set  of  reverse  edges  of  Q  obtained  by 
reversing  the  order  of  all  the  pairs  in  £.  The  reverse  digraph  of  Q ,  denoted  Q ,  is  (U,£,A),  where 
A  =  AT .  The  mirror  digraph  of  Q ,  denoted  Q,  is  (U,  £,  A),  where  £  =  £  U  £  and  A  =  (A  +  A1)/ 2. 
Note  that  L(G)  =  Doni{Q)  -  A{G)  =  DhflG)  -  A(G)T . 

Given  a  digraph  Q  of  order  n,  the  disagreement  function  d>g  :  Mn  — >  M  is  defined  by 


Ui(G)  =  { 


^  aik , 
k=l,k^i 

—dj 


*13) 


®g{x) 


1  X  A 

-  aij(x3  -  xif- 

i,j= 1 


(1) 


To  the  best  of  the  authors’  knowledge,  the  following  is  a  novel  result. 

Proposition  2.2  (Matrix  representation  of  disagreement).  Given  a  digraph  Q  of  order  n,  the 
disagreement  function  :  Mn  — >  M  is  the  quadratic  form  associated  with  the  symmetric  positive- 
semidefinite  matrix 

P(G)  =  \(Dont(G)  +  Din(G)  -  A(G)  -  A(Gf). 

Moreover,  P(G)  is  the  graph  Laplacian  of  the  mirror  graph  G,  that  is,  P(G)  =  L(G)  =  \(^L(G)  + 

m). 
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Proof.  For  x  £  Rn,  we  compute 

x1  P{G)x  =  -xT(Dou t  +  An  -  A  -  A1  )x 

^  n  n  n 

=  -  (  '^2  aijxi  +  ClijXj  —  2  ClijXiXj ) 

*J=1  *J=1  *,J=1 

1  n 

=  —  (  aij(xi  A  Xj  ~  2xixj)) 

i,j=  1 

1  ” 

=  9  A  (A  -  ®t)2  =  $gO). 
v=i 

Clearly  P  is  symmetric.  Since  4>g(x)  >  0  for  all  x  £  Rn,  we  deduce  P(G)  is  positive  semidefinite. 
Since 

n  n  i 

(D(C))ii  =  Q-ij  =  7j(aij  d"  0,ji), 

3=1  3= 1 

we  have  D{G)  =  \(Dont(G)  +  An(6))-  Hence, 

T(S)  =  £>(£)- A(£) 

=  ^(Pout^)  +  An(£0)  -  ^(0)  +^(af)  =  P(0). 

The  last  inequality  follows  from  the  definitions  of  reverse  and  mirror  graphs. 

Remark  2.3.  Note  that  in  general,  P{G)  A  L(Q).  However,  if  the  digraph  Q  is  balanced,  then 
Aut(<?)  =  Dm(Q),  and  therefore, 

$g(x)  =  \xT(D0ut(G)  +  Din(G))x  -  ^ xT(A(G )  +  A(G)t)x 
=  xT  Dout(G)x  —  xT  Ax  =  xTL(G)x- 

This  is  the  result  usually  presented  in  the  literature  on  undirected  graphs.  • 

2.2  Averaging  plus  connectivity  achieves  consensus 

To  each  node  i  £  U  of  a  digraph  Q,  we  associate  a  state  xt  £  R,  that  obeys  a  first-order  dynamics 
of  the  form 


ii  =Ui,  i  £  {1, . . .  ,n}. 

We  say  that  the  nodes  of  a  network  have  reached  a  consensus  if  Xi  =  xj  for  all  i,  j  £  {1, . . .  ,n}. 
Our  objective  is  to  design  control  laws  u  that  guarantee  that  consensus  is  achieved  starting  from 
any  initial  condition,  while  vn  depends  only  on  the  state  of  the  node  i  and  of  its  neighbors  in  Q,  for 
i  £  {1, . . .  ,n}.  In  other  words,  the  closed- loop  system  asymptotically  achieves  consensus  if,  for  any 
xq  £  Rn,  one  has  that  x(t)  — >  {a(l, . . . ,  1)  |  a  £  R}  when  t  —>  +oo.  If  the  value  a  is  the  average  of 
the  initial  state  of  the  n  nodes,  then  we  say  the  nodes  have  reached  average-consensus. 

We  refer  to  the  following  linear  control  law,  often  used  in  the  literature  on  consensus  (e.g., 
see  [7,  10,  16]),  as  the  averaging  protocol: 


U{  —  ^  ]  djj  ( X j  xif 

3= 1 


(2) 
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With  this  control  law,  the  closed-loop  system  is 


x(t)  =  —L(G)x(t).  (3) 

Next,  we  consider  a  family  of  digraphs  {Gi,  ■  ■  ■  ,Gm }  with  the  same  vertex  set  {1, . . .  ,n}.  A 
switching  signal  is  a  map  a  :  M+  x  ]Rn  — >  {l,...,m}.  Given  these  objects,  we  can  define  the 
following  switched  dynamical  system 

x(t)  =  -L(Gk)x(t), 
k  =  a(t,  x(t)). 

Note  that  the  notion  of  solution  for  this  system  might  not  be  well-defined  for  arbitrary  switching 
signals.  The  properties  of  the  linear  system  (3)  and  the  system  (4)  under  time-dependent  switching 
signals  have  been  investigated  in  [2,  5,  8,  23].  Here,  we  review  some  of  these  properties  in  the 
following  two  statements. 

Theorem  2.4  (Averaging  over  a  digraph).  Let  Q  be  a  digraph.  The  following  statements  hold: 

(i)  System  (3)  asymptotically  achieves  consensus  if  and  only  if  G  has  a  globally  reachable  node; 

(ii)  If  G  is  strongly  connected,  then  system  (3)  asymptotically  achieves  average- consensus  if  and 
only  if  G  is  balanced. 

Statement  (i)  is  proved  in  [23,  Section  2],  Statement  (ii)  is  proved  in  [2,  Section  VII]. 

Theorem  2.5  (Averaging  over  switching  digraphs).  Let  {G i, . . . ,  Gm}  be  a  family  of  digraphs  with 
the  same  vertex  set  {1 , ...  ,n},  and  let  a  :  M+  — >  {1, . . . ,  m)  be  a  piecewise  constant  function.  The 
following  statements  hold: 

(i)  System  (4)  asymptotically  achieves  consensus  if  there  exist  infinitely  many  consecutive  uni¬ 
formly  bounded  time  intervals  such  that  the  union  of  the  switching  graphs  across  each  interval 
has  a  globally  reachable  node; 

(ii)  If  each  Gi,  i  €  {1, . . . ,  m},  is  strongly  connected  and  balanced,  then  for  any  arbitrary  piece- 
wise  constant  function  a ,  the  system  (4)  globally  asymptotically  solves  the  average- consensus 
problem. 

Statement  (i)  is  proved  in  [5,  Section  III  B],  Statement  (ii)  is  proved  in  [2,  Section  IX]. 

2.3  Averaging  protocol  over  a  fixed  acyclic  digraph 

Here  we  characterize  the  convergence  properties  of  the  averaging  protocol  in  equation  (3)  under 
different  connectivity  properties  than  the  ones  stated  in  Theorem  2.4,  namely  assuming  that  the 
given  digraph  is  acyclic. 

We  start  by  reviewing  some  basic  properties  of  acyclic  digraphs.  Given  an  acyclic  digraph  G, 
every  vertex  of  in-degree  0  is  named  source,  and  every  vertex  of  out-degree  0  is  named  sink.  Every 
acyclic  digraph  has  at  least  one  source  and  at  least  one  sink.  (Recall  that  sources  and  sinks  can  be 
identified  by  following  any  directed  path  on  the  digraph.)  Given  an  acyclic  digraph  Q,  we  associate 
a  nonnegative  number  to  each  vertex,  called  depth ,  in  the  following  way.  First,  we  define  the  depth 
of  the  sinks  of  Q  to  be  0.  Next,  we  consider  the  acyclic  digraph  that  results  from  erasing  the  0-depth 
vertices  from  Q  and  the  in-edges  towards  them;  the  depth  of  the  sinks  of  this  new  acyclic  digraph 
are  defined  to  be  1.  The  higher  depth  vertices  are  defined  recursively.  This  process  is  well-posed  as 
any  acyclic  digraph  has  at  least  one  sink.  The  depth  of  the  digraph  is  the  maximum  depth  of  its 
vertices.  For  n,  d  €  N,  let  Snj  be  the  set  of  acyclic  digraphs  with  vertex  set  {1, . . . ,  n}  and  depth  d. 

Next,  it  is  convenient  to  relabel  the  n  vertices  of  the  acyclic  digraph  Q  with  depth  d  in  the 
following  way:  (1)  label  the  sinks  from  1  to  no,  where  no  is  the  number  of  sinks;  (2)  label  the 
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vertices  of  depth  k  from  Ej=o  ni  +  1  to  Ej= o  nj  +  nk,  where  rik  is  the  number  of  vertices  of  depth 
k,  for  k  G  {1, . . .  ,  d}.  Note  that  vertices  with  the  same  depth  may  be  labeled  in  arbitrary  order. 
With  this  labeling,  the  adjacency  matrix  A(Q)  is  lower-diagonal  with  vanishing  diagonal  entries, 
and  the  Laplacian  L(Q)  takes  the  form 


or,  alternatively, 


m 


'  0  0 

-021  E]=!  02 j 

_  oni  an 2 


0 

0 


En—l 

j= 1 


dnj  J 


m 


OnoXno  Or)ox(n— no) 

L21  L22 


(5) 


where  0 kxh  is  the  k  x  h  matrix  with  vanishing  entries,  L21  G  flj(n~n°ixno  ancj  £22  g  j^n-n0xn-n0_ 
Clearly,  all  eigenvalues  of  L  are  non- negative  and  the  zero  eigenvalues  are  simple,  as  their  corre¬ 
sponding  Jordan  blocks  are  lxl  matrices. 


Proposition  2.6  (Averaging  over  an  acyclic  digraph).  Let  Q  be  an  acyclic  digraph  of  order  n  with 
no  sinks,  assume  its  vertices  are  labeled  according  to  their  depth,  and  consider  the  dynamical  system 
x(t)  =  — L{Q)x(t )  defined  in  (3).  The  following  statements  hold: 

(i)  The  equilibrium  set  of  (3)  is  the  vector  subspace 

ker  L{Q)  =  {(a;s,xe)  G  Rn°  x  R"_n°  |  xe  =  —  L2ixs }. 


(ii)  Each  trajectory  x  :  R. 
recursively  by 


Rn  of  (3)  exponentially  converges  to  the  equilibrium  x*  defined 


xj  =  < 


xfi  0), 

Ei—1  * 

j= 1  aijxj 

Ei—1  ■ 

j= 1  aL 


i  G  {l,...,n0}, 
i  G  {ra0  +  1, . . .  ,n}. 


(in)  If  Q  has  unit  depth,  then  the  disagreement  function  is  monotonically  decreasing  along  any 
trajectory  of  (3). 

Proof.  Statement  (i)  is  obvious.  Statement  (ii)  follows  from  the  fact  that  —L22  is  Hurwitz  and  from 
the  equilibrium  equality 


i—  1  2  —  1  2—1 

0  =  aP  ( X*3  -  Xi  )  =  aPX*j  -  (  ap)  ■ 

j= 1  3= 1  3= 1 


Regarding  statement  (iii),  when  the  depth  of  Q  is  1,  the  adjacency  matrix  and  the  out-degree  matrix 
are  equal  to,  respectively, 


Ono  x  no 

O710  x  (22—220) 

Ono  x  220 

Ono  x  (22—220) 

—L21 

0(n— 220 )  X  (22— 22q)_ 

5 

_0(n— 22o)  X?2o 

L22 

where  L21  and  L22  are  defined  in  (5).  Therefore,  we  compute 

Ln  Lj 1 

0(n— ng)xno  0(n— no)  X  (n— no). 


m 


where  in  G  jj«oxn0.  According  to  Proposition  2.2,  we  have 

p(Q)=l2(L{G)  +  L{G)) 
6 


The  evolution  of  4>g  along  a  trajectory  of  x  :  M+  — >  1"  of  (3)  is  given  by 

)))  =  -x(t)T(L(g)Tp(g)  +  p(g)L(g))x(t) 

=  -x(t)T L{g)T L(g)x(t)  -  x(t)T L(g)T L(g)x{t) 

=  -x(t)T L(g)T L(g)x(t)  <  o, 

where  in  the  last  equality  we  have  used  the  fact  that  L{Q)T  L{G)  =  L(G)T L(Q)  =  0nxn-  Note  that 
is  strictly  decreasing  unless  x(t)  G  kev  L(G),  i.e. ,  the  trajectory  reaches  an  equilibrium. 

Remarks  2.7.  ( i )  If  the  digraph  has  a  single  sink,  then  the  convergence  statement  in  part  (ii)  of 

Proposition  2.6  is  equivalent  to  part  (i)  of  Theorem  2-4- 

(ii)  The  block  decomposition  of  L(G)  holds  only  for  digraphs  with  depth  1.  Indeed,  statement  (Hi) 
is  not  true  for  digraphs  with  depth  larger  than  1.  The  digraph  in  Figure  1  is  a  counterexample. 


Figure  1:  For  this  digraph  of  depth  2,  the  Lie  derivative  of  the  disagreement  function  (1)  along  the 
averaging  flow  (3)  is  indefinite. 


2.4  Averaging  protocol  over  switching  acyclic  digraphs 

Given  a  family  of  digraphs  T  =  {Gi, . . . ,  Gm }  with  vertex  set  {1, . . . ,  n},  the  minimal  disagreement 
function  4>p  :  Mn  — >  R  is  defined  by 


4>r(x)  =  min  4>gfc(x)- 


(6) 


Let  I(x)  =  argmin{<t>gfc (x)  \  k  G  {1  we  consider  state-dependent  switching  signals  a  : 

Rn  — >  {1, . . .  ,m}  with  the  property  that  a(x)  G  I(x),  that  is,  at  each  x  G  Rn,  a(x)  corresponds  to 
the  index  of  a  graph  with  minimal  disagreement.  Clearly,  for  any  such  a,  one  has  <J>p(x)  =  <3?^^  (x). 
Before  giving  our  result,  we  first  point  out  a  helpful  fact. 

Lemma  2.8.  Let  V  =  {<7i, . . . ,  Gmj  C  <Sn,i.  If  Uke{1,...,m}  Gk  €  Sn,i,  then  for  any  i,j  e  {1, ...  ,m}, 
we  have 

U9i)TL{Qj)  =  0  nxn- 

Proof.  Let  G  =  _  Gk ■  Since  G  G  <Sn,i,  so  we  have,  by  proper  ordering  of  the  nodes, 


L(G) 


OnoXno  Onox(n— no) 

Lr2 1  L22 


L(G) 


L^(n— no)  xno 


L 


T 

21 


0 


(n— no)x(n— no)J 


For  any  i  G  {l,...,m},  Gt  is  a  subgraph  of  G,  so  that  L[Gf)  and  L(Gi)  share  the  same  block 
decompositions  as  stated  in  the  last  equation.  Hence,  the  statement  follows  immediately. 
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Proposition  2.9  (Averaging  over  acyclic  digraphs).  Let  T  =  {Q\, . . . ,  Qm}  C  and  assume  that 
Gk  £  <Sn,  1  and  that  a  :  R™  — >  {1, . . . ,  m }  satisfies  a(x)  £  I(x).  Consider  the  discontinuous 
dynamical  system 

x(t)  =  Y(x(t))  =  ~ L{Gk)x(t ),  for  k  =  a(x(t)),  (7) 

whose  solutions  are  understood  in  the  Filippov  sense.  The  following  statements  hold: 

(i)  The  point  x*  £  Rn  is  an  equilibrium  for  (7)  if  and  only  if  for  each  i  £  I(x*),  there  exists 
scalars  A,;  >  0  and  Ylieit x*)  A*  =  1,  such  that 

x*  £  ker  ^  Y  A .*£(&)) . 

iei(x*) 

In  particular,  if  I(x*)  contains  only  one  element  k*  £  {1  then  (8)  can  be  simplified 

to 


x*  £  ker  L{Qk*).  (9) 

(ii)  Each  trajectory  x  :  R+  — >  Rn  of  (7)  converges  to  the  set  of  equilibria. 

(Hi)  The  minimum  disagreement  function  Tp  is  monotonically  non-increasing  along  any  trajectory 
x  :  R+  —*  Rn  of  (7). 

Proof.  We  investigate  first  smoothness  of  d>p.  Because  — <hp  is  the  maximum  of  the  smooth  functions 
— &gk,  by  Proposition  A. 3,  we  know  that  Tp  is  locally  Lipschitz  and  has  generalized  gradient 

<9<hr(x)  =  co{2 P(Gi)x  j  i  £  I(x)}. 

Let  a  £  £y<h r(x),  then  by  definition,  there  exists  u  =  ~Yliei(x a  xip{Gi)x,  where,  for  each  i  £ 
I(x),  A i  >  0  and  Yliei(x)  A*  =  1>  such  that  a  =  lot(  for  all  (  £  3d>r(.T).  In  particular,  for 
C  =  J2iei(x)  2A iP{Qi)x  £  d$r(x),  we  have 

a=("  E  ^MQi)x)T  [  Y  2XiP&i)x) 

i£l(x)  i£l(x) 

=  -xT(  Y  xMQi))T (  Y  xi(HQi)  +  L(Qi)))x 

i£l(x)  i£l(x) 

=  -xT(  Y  xMQi))T {  Y  xMGi))x-xT(  Y  xMGi)T)(  Y  xMGi))x 

iEl(x)  i£l(x)  i£l(x)  i£l(x) 

=  -xT(  Y  xMGi))T (  Y  xMGi))x<0, 

i£l(x)  i£l(x) 

where  in  the  last  equality  we  have  used  Lemma  2.8.  Moreover, 

a  =  0  x  £  ker  ^  Y^  A iL{Qf)^. 

iei(x) 

In  particular,  if  x  is  not  at  any  switching  surface,  then  I (x)  is  a  set  with  only  one  element  k  £ 
{1, . . . ,  m}  and  chl>r(x)  =  2 P(Qk)x.  Therefore,  Cy$t{x)  =  0  if  and  only  if  x  £  ker  L(Qk ).  Therefore, 
we  conclude  that  for  x  £  Rn  and  a  £  £y3>p(a;),  we  have  a  <  0,  i.e. ,  max£y<I>r(x)  <  0  and 
statement  (i)  is  true.  Resorting  to  the  LaSalle  Invariance  Principle  (Theorem  B.2),  we  deduce  that 
any  trajectory  x  :  R+  — >  R”  of  (7)  converges  to  the  set  of  equilibria  as  stated  in  statement  (i)  and 
statement  (iii)  is  clear. 


Remarks  2.10.  •  Statement  (ii)  in  this  theorem  is  weaker  than  statement  (ii)  in  previous  one 

in  three  ways:  first,  we  are  not  able  to  characterize  the  limit  point  as  a  function  of  the  initial 
state.  Second,  we  require  the  depth  1  assumption,  which  is  sufficient  to  ensure  convergence, 
but  possibly  not  necessary.  Third,  we  establish  only  convergence  to  a  set,  rather  than  an 
individual  point.  It  remains  an  open  question  to  obtain  necessary  and  sufficient  conditions  for 
convergence  to  a  point. 

•  Although  the  statement  (ii)  is  obtained  only  for  digraphs  of  unit  depth,  this  class  of  graphs  is 
of  interest  in  the  forthcoming  sections.  • 


3  Discrete  coverage  control 

In  this  section,  we  first  review  the  multi-center  optimization  problem  and  the  corresponding  coverage 
control  algorithm  proposed  in  [18].  We  then  study  the  multi-center  optimization  problem  in  discrete 
space  and  derive  a  discrete  coverage  control  law.  This  leads  to  a  geometric  object  called  the  discrete 
Voronoi  graph.  Finally,  we  show  that  the  discrete  coverage  control  law  is  an  averaging  algorithm 
over  a  certain  set  of  acyclic  digraphs.  Discrete  locational  optimization  problems  are  discussed  in 
[24,  25,  26], 

We  will  consider  motion  coordination  problems  for  a  group  of  robots  described  by  first  order 
integrators.  In  other  words,  we  assume  that  n  robotic  agents  are  placed  at  locations  pi, . . .  ,pn  G  M2 
and  that  they  move  according  to 


Pi  =  Ui,  *  G  {1, . . . ,  n}.  (10) 

We  denote  by  P  the  vector  of  positions  (p\ ,.. .  ,pn)  G  (M2)n.  Additionally,  we  define 

Scoinc  =  {(pi>  ■  ■  ■  ,Pn)  G  (M2)n  I  Pi  =  Pj  for  some  i  j}, 

and,  for  P  ^  <Scojnc,  we  let  denote  the  Voronoi  partition  generated  by  P,  we  illustrate 

this  notion  in  Figure  2  and  refer  to  [24]  for  a  comprehensive  treatment  on  Voronoi  partitions. 


Figure  2:  The  Voronoi  partition  of  a  rectangle  in  the  plane.  We  depict  the  generators  p\, . . .  ,pn 
elevated  from  the  plane  for  intuition’s  sake. 


3.1  Continuous  and  discrete  multi-center  functions 

In  this  section  we  present  a  class  of  locational  optimization  problems  in  both  continuous  and  dis¬ 
crete  settings.  It  would  be  possible  to  provide  a  unified  treatment  using  generalized  functions  and 
distributions,  but  we  avoid  it  here  for  concreteness’  sake. 

Let  Q  be  a  convex  polygon  in  M2  including  its  interior  and  let  </>  :  M2  — >  be  a  bounded 

and  measurable  function  whose  support  is  Q.  Analogously,  let  {q\, . . . ,  gw}  C  K2  be  a  pointset 
and  {</>i,...,  }  be  nonnegative  weights  associated  to  them.  Given  a  non-increasing  function 

/  :  M+  — *  K,  we  consider  the  continuous  and  discrete  multi-center  functions  TL  :  (M2)n  — ►  M  and 
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Wdsc.rt  :  (R2)n  — >  R  defined  by 


TL{P)  =  [  max  f(\\q  -  Pi\\)(f>(q)dq, 

N 

Wdscrtl-P)  =  PaX  ,  fafiUj  ~Pi  ID- 

J  -  '  U . "> 

Now  we  define 

Sequid  =  {(Pl>  •  •  •  ,Pn)  €  O^2)"-  I  h  ~  Pi\\  =  \W  ~  Pk\\  =  d(q ) 

for  some  q  £  {(71, . . . ,  qjy}  and  for  some  i  /  k}, 

where  d(q)  =  minjer1;  ni.  \\q  —  pj\\.  In  other  words,  if  P  fL  e>Cquid>  then  no  point  q.j  is  equidistant 
to  two  or  more  nearest  robots.  Note  that  5equid  is  a  set  of  measure  zero  because  it  is  the  union  of 
the  solutions  of  a  finite  number  of  algebraic  equations.  Using  Voronoi  partitions,  for  P  0  <Scoinc,  we 
may  write 


=  f(h-Pi\\)(t>{q)dq, 

1^1  JVi(P) 


^dscrt  (P)  ~ 


i=  1  qjeVi(P) 


card (qj,  P ) 


fihj-PiW) 


<i>jf(H-Pi\\)+  J2 


i= 1 


'q^Vi(P) 


qjddVi(P) 


card((/j,  P ) 


fihj-PiW) 


where  card  :  R2  X  (R2)n  — ►  {1, . . .  ,n}  denotes  the  number  of  indices  k  for  which  || q.j  —  p^ ||  = 
minjg{1)  nj  || qj  —  Pi\\.  It  is  easy  to  see  that  card  is  distributed  over  the  Voronoi  graph  and  if  qj  is 
a  point  in  the  interior  of  Vi(P),  then  card (qj,P)  =  1.  For  P  0  5coinc  U  5equ;d,  we  have 


n 

Tfdscrt  (P)  =  fijfiWlj  ~  Pi\\)- 

i= 1  o 

fc-eVi(P) 


Remark  3.1.  T/ie  function  f  plays  the  role  of  a  performance  function.  If  {pi, . . .  ,pn}  are  the 
locations  of  n  sensors,  and  if  events  take  place  inside  the  environment  Q  with  likelihood  <f,  then 
f(\\q  —  Pi||)  is  the  quality  of  service  provided  by  sensor  i.  It  will  therefore  be  of  interest  to  find 
local  maxima  for  Ti  and  Wdscrt-  These  types  of  optimal  sensor  placement  spatial  resource  allocation 
problems  are  the  subject  of  a  discipline  called  locational  optimization  [24,  25,  18].  • 

The  following  result  is  discussed  in[19]  for  the  continuous  multi-center  function. 

Proposition  3.2  (Partial  derivatives  of  H).  If  f  is  locally  Lipschitz,  then  Ti  is  locally  Lipschitz  on 
Qn.  Further,  if  f  is  differentiable,  then  H  is  differentiable  on  Qn\Sco-mc,  and,  for  each  i  £  {1, . . . ,  n}, 

jr(P)  =  f  Jr-fiWQ-pMWdq- 

dpi  JVi(P)  dpi 

We  obtain  the  corresponding  properties  of  7-fdscrt  via  nonsmooth  analysis  as  the  following  propo¬ 
sition. 
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Proposition  3.3  (Generalized  gradient  of  Tfdscrt)-  If  f  is  locally  Lipschitz,  then  "Hdscrt  is  locally 
Lipschitz  on  Qn.  Further,  if  f  is  differentiable,  then  Hdscrt  is  regular  on  Qn  and  its  generalized 
gradient  satisfies 

n  3 

dndscrt(P)  =  E^co{ap/dl®  -Pfcll)  I  k  G  Hoj,P)}> 

3= 1 

where  I(qj,P)  is  the  set  of  indices  k  for  which  f(\\qj  —  Pk ||)  =  mStXie{i,...,n}  f(\\Qj  ~  Pi II);  an<I 
particular,  if  P  ^  5COmC  U  <Scquid,  then  Tidscvt  is  differentiable  at  P,  and  for  each  i  G  {1, . . . ,  n} 

-Qp.Hdscrt{P)  =  22  f>3f^p.f(\\(ij~Pi\\)- 

qj&ViiP) 


Proof.  We  re-write  here  Lfdscrt  as 


N 

Hdscrt  (P)  =  V  max  4>jf(\\qj 

— '  iG{l,...,n} 

3= 1 


N 

Pill  )  =  E^(p)> 

i=i 


where,  for  each  j  G  IV}, 


*j(-P)  =  .  ™ax  Jihj  ~Pi ID- 


We  first  investigate  the  smoothness  of  Fj(P).  By  Proposition  A. 3,  it  is  easy  to  see  that  if  /  is  locally 
Lipschitz,  then  Fj(P)  is  locally  Lipschitz  on  Qn,  so  is  IIdscrt{P)- 

Additionally,  if  /  is  differentiable,  then  Fj(P )  is  regular  on  Qn ,  with  generalized  derivative 

r\ 

3Fj(P)  =  co{—f{\\qj  -pi||)  |  i  G  /(<?j,P)}, 

where  I(qj,P)  is  the  set  of  indexes  k  for  which  f(\\qj  —  Pk ||)  =  Pj{P)-  Since  Hdscrt(-F’)  is  a  finite 
sum  of  Fj(P)  with  nonnegative  weights  (j)j,  so  T^dscrt (P)  is  also  regular  (cf.  Proposition  A. 2)  on  Qn, 
with  the  regularity  of  Fj(P),  we  obtain  further  the  generalized  gradient  of  Hdscrt (P)  as 

N  Q 

dHdscrt (P)  =  22^co{gp^  11^'  “PfclD  I  k  6 

3= 1 


The  expression  for  the  partial  derivative  away  from  5COmc  U5equid  is  easy  to  see. 

Let  diHdscrt(P)  denote  the  ith  block  component  of  <9Hdscrt (P),  the  following  result  is  a  conse¬ 
quence  of  Proposition  3.3. 

Corollary  3.4.  If  f  is  differentiable,  then  for  each  i  G  {1, . . .  ,n}, 


diHdscrt(P)C  22  <f>j7rf(\\qj-Pi\\)  + 

o  °Pi 

qjeVi(P) 


Z ^  ' 

qjedVi(P) 


l 


U 


OPi 


For  particular  choices  of  /,  the  multi-center  functions  and  their  partial  derivatives  may  simplify. 
For  example,  if  f(x )  =  —  x2,  the  partial  derivative  of  the  multi-center  function  Tt  reads  (for  P  ^ 


—  (P)  -  2 My.(P)(Cv.(p)  -  pi), 
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where  mass  and  the  centroid  of  W  C  Q  are 


1 


MW  =  I  4>(q)  dq,  Cw  =  I  qfi(q)dq. 

Jw  Mw  Jw 

Additionally,  the  critical  points  P*  of  Ti  have  the  property  that  p*  =  Cy^p*),  for  i  G  {1, . . . ,  n};  these 
are  called  centroidal  Voronoi  configurations.  Analogously,  if  f(x)  =  —x2,  the  discrete  multi-center 
function  Hdscrt  reads 

N 

Hdscrt (P)  =  “XX  ■  ?laX  -  Pill2; 


and  its  generalized  gradient  is 


N  Q 

dHdscrt(P)  =  XX 4>]  co |2(^  -Pk)~^p  |  ^  e  i(qj,P) }• 


i=i 

For  each  j  €  {1, . . . ,  AT},  assume  the  scalars  A y ,  i  G  I(qj,  P),  satisfy 

Ay  ’  0,  ^  (  Ay  —  1. 

i£l(qj,P) 

Next,  define  (Mdscrt)yi(p)  and  (Cdscrt)y(p),  respectively,  as 

(Mdscrt)yj(P)  =  XX  $3  XX  At?<Ai  =  X]  Ab<Ai, 

®eVi(P) 


(ID 


(12) 


v/vo 

qjedVi(P)  qj£Vi(P) 


(Cdscrt)y(P)  —  \ 


Pi, 


1 


(A^dscrt)  Vi(P) 


XX  Ad 

«,-eVi(P) 


if  (Mdscrt)Vi(P)  —  0, 
otherwise. 


dP’ 


Lemma  3.5.  Given  f{x)  =  — x2,  P*  is  a  critical  point  of  d Hdscrt,  *.e.,  0  G  dHdSCrt(P*)i  if  kind 
only  if  for  any  j  G  {1,...,1V},  there  exist  Xij  as  in  (12),  such  that  p*  =  (CdSCrt)y(p*),  /or  eac/i 
i  G  {1, . . .  ,n}. 

Proof.  Given  scalar  satisfying  (12),  define 

N 

™  =  XX  ^  XX  2Ao  (h  -  pD  1 

i=l  kei(qj,P) 

then  it  is  clear  that  u;  G  <9Hdscrt(P*)-  Let  Wi  denotes  the  ith  component  of  w,  since 

Wi  =  2  X:  ~Pi)  =  ^  XX  ~~  ( Cdscrt ) V< (P* ) )  =  0, 

«ev4(p*)  «aeVi(p*) 

so  rc  =  0  and  0  G  <9iFdscrt(P*)- 

On  the  other  hand,  if  P*  is  a  critical  point,  then  there  exists  scalars  Ay  satisfying  (12),  such 
that 

w  =  XX  h  XX  2A0 (Qj  ~  Pk)  =  °> 

J=1  keI(qj,P) 

which  implies,  for  each  i  G  {1, . . . ,  n} 

uh  =  2  XX  Ai./°iL/,  -P*)  =  0, 

qj£Vi(P*) 

Solve  this  linear  equation,  we  obtain 

Pi  =  ( Cdscrt)y(P*) ,  i  €  {1,  •  •  ■  , 
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3.2  Continuous  and  discrete  coverage  control 

Based  on  the  expressions  obtained  in  the  previous  subsection,  it  is  possible  to  design  motion  coor¬ 
dination  algorithms  for  the  robots  p\, ...  ,pn.  We  call  continuous  and  discrete  coverage  control  the 
problem  maximizing  the  multi-center  functions  TL  and  Ltdscrt;  respectively.  The  continuous  prob¬ 
lem  is  studied  in  [18].  We  simply  impose  that  the  locations  pi, . . .  ,pn  follow  a  gradient  ascent  law 
defined  over  the  set  Qn  \  Sco inc-  Formally,  we  set 

Ui  =  kprop^iP),  (13) 

opi 

where  kprop  is  a  positive  gain.  Note  that  this  law  is  distributed  in  the  sense  that  each  robot  only 
needs  information  about  its  Voronoi  cell  in  order  to  compute  its  control. 

For  discrete  coverage  control,  we  adopt  the  following  discontinuous  control  law,  for  each  robot 
i  G  {1, . .  .,n} 


—  1'prop  XjfPf 


where  X%  :  Qn  — ►  M2  is  defined  as 

*i(p)=  E 

QjeVi(P) 


h  9 

card (qj,  P )  dp , 


(14) 


Note  that  X%  is  continuous  at  P  G  Qn  \  <Scoinc  U5equid,  and  satisfies 


Xi(P) 


d74dscrt 

dpi 


(P)- 


Like  control  law  (13),  the  discontinuous  control  law  (14)  is  also  distributed.  Define  the  vector  field 
X  =  [X\,  X-2, . . . ,  Xn]7 ,  we  have 

P  =  kpmpX(P).  (15) 

Since  X(P)  is  discontinuous  at  P  €  5coinc  U  5equid,  we  understand  the  solution  of  this  equation  in 
the  Filippov  sense  following  [27],  and  the  existence  of  Filippov  solution  is  guaranteed.  We  then 
investigate  the  properties  of  the  solution  and  analysis  the  convergence  of  (13)  and  (14). 

Proposition  3.6  (Continuous  coverage  control;  [18,  19]).  For  the  closed-loop  systems  induced  by 
equation  (13)  starting  at  Pq  €  Qn  \  Sco\nc,  the  agents  location  converges  asymptotically  to  the  set  of 
critical  points  of  PL. 


Proposition  3.7  (Discrete  coverage  control).  For  the  closed-loop  systems  induced  by  equation  (14) 
starting  at  Pq  €  Qn  \  Sco inc,  the  agents  location  converges  asymptotically  to  the  set  of  critical  points 
°f  ^dscrt  • 


Proof.  Note  that 

K[kpropX](P)  = 

^prop^Ttdscrt  ( P )  • 

Given  this  property,  the  following  proof  is  essentially  the  same  as  the  proof  of  Proposition  2.9  in  [28]. 
We  refer  the  interested  reader  to  [28]  for  technical  details. 


3.3  Discretizing  continuous  settings 

In  this  section  we  discuss  the  relationship  between  the  discretization  of  continuous  locational  opti¬ 
mization  problems  and  discrete  locational  optimization  problems. 

As  before,  let  Q  be  a  convex  polygon  in  M2  including  its  interior,  and  let  <f>  :  M2  — >  M+  be  a 
bounded  and  measurable  function  whose  support  is  Q.  We  shall  consider  a  sequence  of  pointsets 
{Qii  •  •  •  j  QNk}keN  C  M2  and  of  nonnegative  weights  {4>\,  ■  ■  ■ ,  }feeN-  Accordingly,  we  can  define  a 


sequence  of  discrete  multi-center  functions  H^scrt,  for  k  E  N.  The  sequence  {q\, . . . ,  q%k}keN  C  R2 
is  dense 1  in  Q  if,  for  all  q  £  Q, 


lirn  min{||g  -  z\\  \  z  E  {q\, , . . ,  qkN  }}  =  0. 
fc^+oo 

Given  a  pointset  qi,...,qjg,  let  V(qi, . . . ,  qjg)  denote  the  Voronoi  partition  it  generates  and 
define  the  associated  weights 

fij  =  [  0{q)dq.  (16) 

Proposition  3.8  (Consistent  discretization).  Assume  that  f  is  continuous  almost  everywhere,  that 
the  sequence  {q\,  ■  ■  ■  ,q%  }keN  C  M2  is  dense  in  Q,  and  that  the  sequence  of  weights  are  defined 
according  to  (16).  Then  {'M^scrt  converges  pointwise  to  H,  that  is,  for  all  P  E  Qn  , 

lim  Hkdscit(P)  =  H(P). 

k— >+oo 

Additionally,  if  f  is  continuously  differentiable,  then  for  P  E  Qn  \  <Scoinc  and  each  i  E  {1 ,n}, 
any  sequence  Xk  E  <9,;?4(}scr1:(-P),  k  E  N,  satisfies 

lim  xk  =  ^{P). 

k^+oo  dpi 

Proof.  For  k  E  N,  given  the  pointset  {q\, . . . ,  (fifi  },  we  define  the  projection  projfc  :  Q  —*■  {q\, . . . ,  q^  } 


projfc(g)  =  argmin{||g  -  z\\  \  z  E  {q*, qNk}}- 
Because  of  the  vanishing  dispersion  property,  we  know  that,  for  all  q  E  Q, 


lim  projfc(g)  =  q. 

«— >+o O 


Therefore,  we  compute 


^dscrt(p)  =  }nax  ,  /(IkJ  -  Pill)  /  <Kv)dq 

Nk 

=  Y]  /  max  .  f(hj  -  Pi\\)4>{q)dq 

~{Jvj(ql...,qkNk)ie{l,-,n} 

=  /  max  f(\\pro}k{q)  ~  Pi\\)(f>{q)dq. 


Because  /  is  continuous  almost  everywhere,  we  have 

lim  nkdscrt(P) 

k — »-|-oo 


=  lim  /  max  /( ||  proj N  (q)  -  Pi\\)(f>(q)dq 
=  /  max  f(||  lim  proj^ (q)  -  Pi\\)fi{q)dq 

jQie{l,...,n}  k^+oo 

=  1  max  f(\\q-Pi\\)(j){q)dq  =  n(P). 

Jq  »6{l,...,n} 

1This  is  equivalent  to  asking  that  the  sequence  has  vanishing  dispersion ;  the  dispersion  of  a  pointset  {gi, . . . ,  f/jv} 
in  the  compact  set  Q  is  max  min  I  In  —  2 II. 

qeQ  ze{qi,...,qN}" 
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Define 


diWdscAP)  =  Y  -^11)  +  <^co{ 


tfcViiP) 


q«£dVi(P) 


By  Corollary  3.4,  if  /  is  differentiable,  then  diH\scrt(P)  C  d^H^SCTt(P) .  Suppose  x*k  G  d*H^scrt(P), 
then  there  exists  scalars  A y  G  [0, 1],  such  that 


xk  =  Y  <f>jjr-f(\\q!j-Pi\\)+  Y  -Pi\ 


dpi' 


q*eVi{P) 


q^edVi(P) 


dpi' 


(18) 


Substitute  (16)  into  (18),  we  obtain 


* 

xk  = 


E 

qHVi(P) 


-  Piling 


+  y[ 

q^edVi(P)  “'E/(<?i.-.9at) 


Aij^>(g)^:/(||gj:  -Pi||)dg- 


Since  /  is  continuously  differentiable,  so  for  P  G  \  5Coinc>  we  have 


lim  xt  = 

k—>-\-oo 


'v 


&{q)J>-f{\\q-Pi\\)dq  + 

i(p)  dPi 


=  |V). 

Opi 


Hence, 


lim  xk  =  ^(P). 

k^+o o  Opi 


3.4  The  relationship  between  discrete  coverage  and  averaging  over  switching 
acyclic  digraphs 

As  above,  let  Q  be  a  convex  polygon,  let  {pi,...,pn}  C  Q  be  the  position  of  n  robots,  let 
{<71, . . . ,  qi\r}  C  Q  be  N  fixed  points  in  Q  with  corresponding  nonnegative  weights  {0i, . . . ,  (/>n}, 
and  let  I(qj,P)  be  the  set  of  indices  k  for  which  \\q3  —  Pk\\  =  minie{i,...,n)  | q.j  —  Pi\\.  We  begin  by 
defining  a  useful  digraph  and  a  useful  set  of  digraphs. 

A  discrete  Voronoi  graph  f/dscrt-Vor  is  a  digraph  with  (n  +  N)  vertices  {pi, . . .  ,pn,  qi, . . . ,  </tv}, 
with  N  directed  edges 

{(pi,  qj)\  for  each  j  G  {1, . . . ,  N},  pick  one  and  only  one  i  G  I(qj,  P)}, 

and  with  corresponding  edge  weights  (f>j ,  for  all  j  G  {1, . . . ,  N}.  We  illustrate  this  graph  in  Figure  3. 
With  our  definition,  it  is  possible  for  one  vertex  set  to  generate  multiple  discrete  Voronoi  graphs. 
We  will  denote  the  nodes  of  f/dscrt-Vor  by  Z  =  (z\, . . . ,  zn+jq)  €E  (M2)n+Ar,  the  weights  by  aap,  for 
a,  /3  G  {1, . . . ,  n  +  N},  with  the  understanding  that: 

_  I  Pa,  if  a  G  {1, . . .  ,n}, 

1  qa-n,  otherwise, 

and  that  the  only  non- vanishing  weights  are  aap  =  4>j  when  /3  =  n  +  j,  for  j  G  {1, ...  ,1V},  and 
when  a  G  {1, . . . ,  n}  corresponds  to  the  robot  pa  closest  to  qj  and  ( pa ,  qj)  is  a  directed  edge  of  the 
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Figure  3:  The  discrete  Voronoi  graph  over  3  robots  and  6x9  grid  points.  This  illustration  is  to  be 
compared  with  the  Voronoi  partition  illustrated  in  Figure  2.  The  edges  have  top/down  direction. 


graph  !?dscrt-Vor ■  Note  that  f/dscrt-Vor  depends  upon  Z.  Since  {(?i, . . . ,  gjv}  C  Q  are  fixed,  when  we 
need  to  emphasize  this  dependence,  we  will  simply  denote  it  as  f/dscrt-Vor(-P)- 

Let  us  now  define  a  set  of  digraphs  of  which  the  discrete  Voronoi  graphs  are  examples.  Let 
F(N,n )  be  the  set  of  functions  from  {l,...,iV}  to  {l,...,n}.  Roughly  speaking,  a  function  in 
F(N,  n)  assigns  to  each  point  qj,  j  E  {1, ... ,  N},  a  robot  pi,  i  E  {1, . . . ,  n}.  Given  h  E  F(N,  n),  let 
Qh  be  the  digraph  with  (n  +  N)  vertices  {pi, . . .  ,pn,  q±, . . . ,  <pv},  with  N  directed  edges 

{(Ph(j)i  Qj)}je{i,...,N}i 

and  corresponding  edge  weights  <j)j,  j  E  {1, . . . ,  N}.  With  these  notations,  it  holds  that  f/dscrt-Vor(T>)  = 
Qh*(;P)  with  any  function  h*  :  {1, . . . ,  N}  x  Qn  — >  {1, . . . ,  n}  which  satisfies 

h*{j,P)  G  argmin{||gj  -  pi\\  \  i  E  {1, . . .  ,n}}. 

Let  us  state  a  useful  observation  about  these  digraphs. 

Lemma  3.9.  The  set  of  digraphs  Qh,  h  E  F(N,n),  is  a  set  of  acyclic  digraphs  with  unit  depth, 
i.e.,  it  is  a  subset  of  Su+n,i  (see  definition  in  Subsection  2.3).  Moreover,  ^h&F(N,n)  Gh  is  an  acyclic 
digraph  with  unit  depth,  i.e.,  U heF(N,n)  Gh  G  <Sn+iV,i- 

For  h  E  F(N,n),  let  us  study  appropriate  disagreement  functions  for  the  digraph  Qh-  We  define 
the  function  Qgh  :  ( M?)n+N  — >  R  by 

1  n+N 

^Qh(Z)\z=(p1,...,pn,q1,...,qN)  =  X  XT  aa/3|lzo  —  z&\\ 

a, 13=1 

1  N 

=  2^2  hH  -PhU)\\2’ 

W  3= 1 

because  the  weights  aap,  a,  (3  E  {1, . . . ,  n  +  N}  of  the  digraph  Qh  all  vanish  except  for  ah(^  j  =  <pj, 
j  E  { 1 , . . . ,  A  } . 

We  are  now  ready  to  state  the  main  result  of  this  section.  The  proof  of  the  following  theorem 
is  based  on  simple  book-keeping  and  is  therefore  omitted. 

Theorem  3.10  (Discrete  coverage  control  and  averaging).  The  following  statements  hold: 

(i)  The  discrete  multi- center  function  Lfdscrt  with  f(x)  =  —  x2 ,  and  the  minimum  disagreement 
function  over  the  set  of  digraphs  Qh,  h  E  F(N,n),  satisfy 

1  1  N 

-xHdscrt(-P)  =  o^2.  ™in  ,  (hhj  ~  Pi  f 

z  2  r-*  ie{i,...,n} 

3=1 

!  N 

=  2*52  fall®  ~Ph*(j)\\2 
3=1 

^  ^dscrt- Vor  (pl,...,pn,q1,...,qN) 

=  min  <f>gh(pi,...,pn,qi,---,qN)- 

heF(N,n ) 
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(ii)  For  P  0  <Scoinc  U  tSequid?  discrete  coverage  control  law  for  f(x)  =  —  x2  and  the  averaging 
protocol  over  the  discrete  Voronoi  digraph  satisfy,  for  i  £  { 1 , . . . ,  n}, 


l<9Hdscrtfp> 

2  dPi  [  ’ 


n+N 

Moj  Pi)  =  a*p(zh  ~  2«), 

qjeVi(P)  (3=1 


where  za  and  aag,  a,  (3  £  {1, . . . ,  n  +  N},  are  nodes  and  weights  of  £/dsCrt- vor  •  Accordingly, 
the  discontinuous  coverage  control  system  (15),  for  f(x)  =  — x2,  and  the  averaging  system  (7) 
over  the  set  of  digraphs  Gh,  h  £  F(N,  n )  satisfy,  for  i  £  {1 , ,n}, 

l-K[Xi]{P)  =  K[Yi]{Z), 

where  Z  =  (pi, . . .  ,pn,  qi, . . . ,  q tv),  X,  and  Y \  are  the  ith  2-dimensional  block  component  of  X 
and  Y ,  respectively. 

(Hi)  Any  P*  £  Qn  is  an  equilibrium  of  the  discrete  coverage  control  system  with  f(x)  =  —x2  if  and 
only  if  Z*  =  (p\, . . .  ,p* ,  q\, . . . ,  q n)  is  an  equilibrium  of  system  (7)  over  the  set  of  digraphs 
Qh,  h  £  F(N,n),  that  is: 


Vj  £  {1, . .  .,N},  3X  ij  as  in  (12),  such  thatp*  =  (Cds  crt)Vi(p*),  Vi  £  {1,» . .  ,n}, 

3/ik  >  0  and  (ik  =  1,  such  that  Z*  £  ker  E  hkL(G dscrt-Vor(-^*)))  > 
k  k 

where  { 0dscrt- Vor  ( -^ * )} k  are  possible  discrete  Voronoi  graphs  generated  by  Z* . 

(iv)  Given  any  initial  position  of  robots  Pq  £  Qn,  the  evolution  of  the  discrete  coverage  control 
system  (15)  and  the  evolution  of  the  averaging  system  (7)  under  the  switching  signal  a  : 
Qn  —■ ►  {Gh  I  h  £  F(N,  71)}  defined  by  cr(P)  =  (5dscrt-Vor(-^)  are  identical  in  the  Filippov  sense 
and,  therefore,  the  two  systems  will  converge  to  the  same  set  of  equilibrium  placement  of  robots, 
as  described  in  (Hi). 


3.5  Numerical  simulations 

To  illustrate  the  performance  of  the  discrete  coverage  law  as  stated  in  Proposition  3.7  and  to 
illustrate  the  accuracy  of  the  discretization  process,  as  analyzed  in  Proposition  3.8,  we  include  some 
simulation  results.  The  algorithms  are  implemented  in  Matlab  as  a  single  centralized  program. 
As  expected,  the  simulations  for  the  discrete  coverage  law  are  computationally  intensive  with  the 
increase  in  the  resolution  of  the  grid.  We  illustrate  the  performance  of  the  closed-loop  systems  in 
Figures  4,  5,  6  and  7. 

4  Conclusions 

We  have  studied  averaging  protocols  over  fixed  and  controlled-switching  acyclic  digraphs,  and  char¬ 
acterized  their  asymptotic  convergence  properties.  We  have  also  discussed  continuous  and  discrete 
multi-center  locational  optimization  functions,  and  distributed  control  laws  that  optimize  them. 
The  main  result  of  the  paper  shows  how  these  two  sets  of  problems  are  intimately  related:  discrete 
coverage  control  laws  are  indeed  averaging  protocols  over  acyclic  digraphs.  As  a  consequence  of  our 
analysis,  it  may  be  argued  that  the  coverage  control  problem  and  the  consensus  problem  are  both 
special  cases  of  a  general  class  of  distributed  optimization  problems. 
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Figure  4:  Continuous  coverage  law  for  32  agents  on  a  convex  polygonal  environment,  with  density 
function  cj)  =  exp(5.(— x2  —  y 2))  centered  about  the  gray  point  in  the  figure.  The  control  gain 
in  (13)  is  kpTop  =  1  for  all  the  vehicles.  The  left  (respectively,  right)  figure  illustrates  the  initial 
(respectively,  final)  locations  and  Voronoi  partition.  The  central  figure  illustrates  the  gradient 
descent  flow.  Figure  taken  from  [18]. 


Figure  5:  Simulation  of  discrete  coverage  law  with  159  grid  points. 


Figure  6:  Simulation  of  discrete  coverage  law  with  622  grid  points. 


Figure  7:  Simulation  of  discrete  coverage  law  with  2465  grid  points. 
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A  Basic  ideas  on  nonsmooth  analysis 

In  this  appendix  we  review  some  properties  in  basic  calculus  of  generalized  gradient  of  functions  [29]. 

Proposition  A.l  (Scalar  multiples).  If  f  :  RM  — >  M  is  locally  Lipschitz  at  x  €  M.m  ,  then,  for  any 
scalar  s,  sf  is  locally  Lipschitz  at  x  and 

d(sf){x)  =  sdf(x). 

Proposition  A. 2  (Finite  sums).  If  fi  :  Mm  — >  R,  i  G  {1, . . .  ,  m},  are  locally  Lipschitz  at  x  £  Mm, 
then,  for  any  scalars  {si, . . . ,  sm},  Y^T=i  sifi  ^oca^U  Lipschitz  and 

m  m 

^2  Sifo)  C  ^2  Si9MX)’ 

i= 1  i= 1 

where  equality  holds  if  each  fi  is  regular  at  x  and  each  Si  is  nonnegative. 
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Proposition  A. 3.  Let  ff  :  M.M  — >  R,  k  G  {1, . . . ,  m}  6e  locally  Lipschitz  functions  at  x  G  and 
let  f(x')  =  max{ fk{x')  \  k  G  {1, . . , ,  m}}.  Then, 

(i)  f  is  locally  Lipschitz  at  x, 

(ii)  if  I(x')  denotes  the  set  of  indexes  k  for  which  fk(x')  =  f(x'),  we  have 

df(x)  C  co {dfi(x)  |  i  G  I(x)}  , 

and  if  fi,  i  6  I{x),  is  regular  at  x,  then  equality  holds  and  f  is  regular  at  x. 

B  Stability  analysis  via  nonsmooth  Lyapunov  functions 

Throughout  the  paper,  we  define  the  solutions  of  differential  equations  with  discontinuous  right- 
hand  sides  in  terms  of  differential  inclusions  [27].  Let  F  :  M.N  — >  2K"Y  be  a  set-valued  map.  Consider 
the  differential  inclusion 


x  G  F(x) .  (19) 

A  solution  to  this  equation  on  an  interval  [to,  ti]  C  M  is  defined  as  an  absolutely  continuous  function 
x  :  [t0,ti]  — * >  such  that  x(t)  G  F(x(t))  for  almost  all  t  G  [to » ^l] -  Given  xo  G  M.N ,  the  existence  of 

at  least  a  solution  with  initial  condition  xo  is  guaranteed  by  the  following  lemma. 

Lemma  B.l.  Let  the  mapping  F  be  upper  semicontinuous  with  nonempty,  compact  and  convex 
values.  Then,  given  xo  G  ,  there  exists  a  local  solution  of  (19)  with  initial  condition  xq. 

Now,  consider  the  differential  equation 

x(t)  =  X(x(t)) ,  (20) 

where  X  :  — >  M.N  is  measurable  and  essentially  locally  bounded.  Here,  we  understand  the 

solution  of  this  equation  in  the  Filippov  sense,  which  we  define  in  the  following.  For  each  x  G  W.N, 
consider  the  set 

A[X](x)  =  f|  p|  co  {X(Bn(x,S)\S)}, 

5>0[i{S)=0 

where  /. i  denotes  the  usual  Lebesgue  measure  in  Mw.  Alternatively,  one  can  show  [30]  that  there 
exists  a  set  Sx  of  measure  zero  such  that 

K[X](x)  =  co{  lim  X(xi)  \  x*  — >  x  ,  x*  fL  S  U  Sx}  , 

2— >+00 

where  S  is  any  set  of  measure  zero.  A  Filippov  solution  of  (20)  on  an  interval  [to,H]  C  1  is 
defined  as  a  solution  of  the  differential  inclusion  x  G  K[X\(x).  Since  the  multivalued  mapping 
K[X\  :  1^  — >  2r  is  upper  semicontinuous  with  nonempty,  compact,  convex  values  and  locally 
bounded  (cf.  [27]),  the  existence  of  Filippov  solutions  of  (20)  is  guaranteed  by  Lemma  B.l. 

Given  a  locally  Lipschitz  function  /  :  — >  M,  the  set-valued  Lie  derivative  of  f  with  respect  to 

X  at  x  is  defined  as 

Cxf{x)  =  {a  G  R  |  G  K[X\(x)  such  that  £  •  v  =  a ,  VC  G  df(x)}  . 

For  each  x  G  M'v,  Cxf(x )  is  a  closed  and  bounded  interval  in  M,  possibly  empty.  If  /  is  continuously 
differentiable  at  x,  then  Cxf{x)  =  {df  ■  v  \  v  G  K[X](x)}.  If,  in  addition,  X  is  continuous  at  x, 
then  Cxf(x)  corresponds  to  the  singleton  {Cxf(x)},  the  usual  Lie  derivative  of  /  in  the  direction 
of  X  at  x. 

The  following  result  is  a  generalization  of  LaSalle  Invariance  Principle  for  differential  equations 
of  the  form  (20)  with  nonsmooth  Lyapunov  functions.  The  formulation  is  taken  from  [31],  and 
slightly  generalizes  the  one  presented  in  [32]. 
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Theorem  B.2  (LaSalle  principle).  Let  f  :  Rw  — >  R  be  a  locally  Lipschitz  and  regular  function. 
Let  xo  G  R^  and  let  /_1(<  /(to),to)  be  the  connected  component  of  {t  G  Rw  |  /( x)  <  /( To)} 
containing  x$.  Assume  the  set  /_1(<  /( To), To)  is  bounded  and  assume  either  ma x£\-/(t)  <  0  or 
Cxf(x)  =  0  for  all  x  G  /_1(<  /( To),  To).  T/ien  /_1(<  /(to),to)  is  strongly  invariant  for  (20).  Let 

ZXJ  =  {x  G  RN  |  0  G  £a7(t)}  • 

Then,  any  solution  x  :  [to,  +oo)  — >  MN  of  (20)  starting  from  to  converges  to  the  largest  weakly 
invariant  set  M  contained  in  Z xjF\f~l(<  /( To),  To).  Furthermore,  if  the  set  M  is  a  finite  collection 
of  points,  then  the  limit  of  all  solutions  starting  at  tq  exists  and  equals  one  of  them. 
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